航天器视场遮挡分析方法研究

敏感器视场遮挡分析是航天器构型布局设计的重要组成部分，针对大型在轨变构型航天器敏感器视场遮挡分析工作量大、出错率高、遮挡率无法量化等特点，本文提出一种基于敏感器第一视角的视场遮挡快速分析方法。该方法通过设置第一视角图像蒙版及图像映射，模拟方锥形、圆锥形、球形等各类敏感器视场，通过逐个像素比对，计算视场遮挡率和遮挡方位。使用此方法建立了空间站敏感器视场遮挡分析系统，结果表明新方法能够显著提高视场遮挡分析效率和准确率，并且为机构类部件运动过程对敏感器视场遮挡的分析提供了有效工具。     

超临界环境下煤油替代物液滴燃烧特性数值研究

为研究液氧/煤油火箭发动机燃烧室内经喷注形成的煤油液滴的燃烧过程，基于实际气体状态方程、高压热物性修正、高压气液平衡和详细化学反应动力学，建立一维的全瞬态液滴燃烧模型，对超临界环境下两组分煤油替代物液滴的燃烧特性及液滴初始直径的影响进行仿真研究。结果表明，在超临界环境下，相比于煤油液滴纯蒸发过程，煤油液滴燃烧过程的迁移时刻大大提前；煤油液滴着火之后很快进入超临界燃烧阶段，此时液滴燃烧过程可以看成中心附近的燃料高浓度区与外侧氧气高浓度区之间的扩散燃烧过程；煤油液滴的火焰半径先增大，达到最大值之后开始减小，并减小为零，火焰温度在着火之后快速上升至最大值，并基本保持不变，在火焰半径减小为零之后开始降低；随着液滴初始直径的增大，火焰特性以及液滴中心参数变化曲线趋势不变、整体延迟，着火时间、迁移时间和液滴寿命增大。     

弹载变转速定量泵高动态压力控制特性分析

传统弹载电液伺服系统的液压能源采用有刷直流电机驱动恒压变量泵的方式，泵源压力在导弹飞行过程中设定值恒定，但在平飞段过程长、负载小，对能源压力需求低，造成弹上能源总利用率不高的问题，提出无刷直流电机驱动小排量液压泵的压力控制泵源的设计，在低气动负载阶段适当降低系统压力，进而降低节流损失与泄漏量，提高能量利用率。分析了变转速定量泵的流量压力特性曲线，通过建立基于变转速压力控制泵源的位置电液伺服系统数学模型与仿真模型，分析了变转速定量泵压力控制系统在典型工况下的位置跟踪、系统压力变化与能量利用效率，验证了压力控制系统能够在导弹平飞段保证位置伺服系统跟踪正常的前提下，能有效提高系统效率，证明了变转速定量泵压力控制系统的有效性与可行性。     

对于动态目标的航天器鲁棒PD+快速跟踪控制方法

摘要：针对卫星对于空间动态目标快速、稳定地跟踪、控制目标，同时考虑平台模型的不确定性、外部随机干扰、系统控制力矩与角速度约束等因素，设计PD+控制器实现对于动态目标的快速、稳定跟踪；在经典PD控制器的基础上设计控制添加项使得系统能够按照既定轨迹运动；采用变结构的手段实现系统收敛速度的提升；合理设计Lyapunov函数的结构，引出角速度、四元数的耦合项对V函数进行改良，简化系统稳定性证明与分析的过程；讨论系统最极端情形，通过对V函数上下界的讨论分析系统该情形下的稳定性；最后通过数值仿真验证所提出算法的有效性与优越性.     

基于Hilbert-Huang变换的航天器发射段力学环境分析

文章基于Hilbert-Huang变换,对某载人航天器发射过程中整流罩分离期间测量的力学环境参数进行处理分析,利用HHT的基函数自适应性特点,得到振动时间历程信号的时频域Hilbert谱,准确识别出信号主频率及对应时间范围。此应用案例可为地面力学试验条件制定以及运载接口条件确定提供参考。     

跨介质航行器流体动力外形仿生设计

跨介质航行器是一种可在空中与水中两栖巡航并能自由穿越水气界面的新概念海空两栖无人运动平台,需要具备良好的入水和出水性能及较小的水下流体阻力。采用两种生物组合仿生的设计思想,利用全自动三维影像扫描仪获取具有优良入水性能的翠鸟头部和龙虱身体的形体特征点云数据,应用傅立叶级数拟合方法拟合特征曲线,设计了一种跨介质航行器流体动力外形,为减小航行器入水冲击载荷、降低水下航行阻力提供了设计方案。     

耦合热、动载荷的超超临界汽轮机迷宫密封动力特性

密封汽流激振严重影响超超临界汽轮机的安全运行,采用DEFINE_CG_MOTION和DEFINE_PROFILE控制宏建立转子的涡动方程,通过Workbench流固耦合方法计算热、动载荷下密封齿形变,根据快速傅里叶变化得到机组运行时的密封动力特性,并对转子稳定性进行分析。结果表明:蒸汽可导致密封齿膨胀变形,温度对密封齿长度变化影响可达1%~1.5%,压力和离心作用对其影响较小。热、动载荷使迷宫密封直接刚度减小,直接阻尼先增加后减小,交叉刚度先减小后增加,动力系数的最大变化为原来的2倍。35~55 Hz内转子稳定裕度急剧下降,转子对密封汽流激振更敏感。热、动载荷引起的压力波动集中在低频范围,密封周向压力波动可增高18.5 kPa。密封高压区的压力波幅剧增是汽流激振显著的主要原因。      

类X-37B升力体再入飞行通信特性中断研究

针对再入非平衡等离子流场的数值仿真方法及等离子体鞘套对不同波段电磁波传输特性的预测开展研究。首先，利用某工程钝头体飞行器再入实测通信飞行测量数据，初步校验数值仿真等离子体特性和通信特性方法的正确性；其次，在典型弹道下针对升力体再入飞行器开展等离子体鞘套数值仿真，分别给出S波段、Ka波段与K波段电磁波衰减特性，研究指出结合当前国内中继星体制情况采用高频段电磁波Ka波段或K波段中继通信方式是解决通信黑障的技术途径之一；最后，研究了电磁波数衰减受传播方向角和攻角的影响规律。研究结果可为发展解决我国再入飞行器通信黑障的设计方案提供参考。     

非预混条件下的旋转爆轰燃烧室双波头演化过程数值模拟

针对旋转爆轰燃烧室双波头演化过程中流场结构变化的问题,对非预混条件下的旋转爆轰燃烧室从起爆到形成稳定的双波头过程进行了数值模拟研究。研究结果表明,从起爆到形成稳定爆轰过程,燃烧室主要经历了起爆、爆轰波对撞和稳定爆轰三个阶段；在爆轰波对撞阶段,首次对撞是两个爆轰波间的对撞,由于对撞点处缺少新鲜混合气,从而在对撞结束后衰减为两个压力波。第二次对撞是两个压力波间的对撞,因为在第二次对撞点附近存在新鲜混合气来支撑爆轰波的持续传播,故对撞结束后产生了一个爆轰波和一个较弱的压力波；第二次对撞发生后,燃烧室内的压力波反射叠加并形成局部高压区,此高压区压缩气体使气体温度升高,高温气体引燃混合气后,最终发展成为第二个爆轰波；稳定阶段,两个爆轰波均能稳定自持传播,爆轰波峰面压力可达1.45MPa,波后温度为2500K,爆轰波速度稳定在1738m/s,产生的推力与比冲分别为79.76N和2312.15s；斜激波的存在使燃烧室出口平面流场产生了较大波动。      

A comprehensive study of the 2016 Mw 6.0 Italy earthquake based on high-rate (10 Hz) GPS data

The 2016 M_w 6.0 Italy earthquake is successfully recorded by the near-field 10?Hz GPS and 200?Hz Strong Motion (SM) stations, providing valuable data for this study. A comprehensive study of this earthquake is carried out based on GPS data, which contains coseismic deformations analysis, noise analysis, seismic wave picking, and magnitude determination. The noise of most GPS-derived displacement waveforms can be described as a combination of white noise, flicker noise, and random walk noise after the earthquake occurrence, and the spectral indices vary significantly for most stations, implying that the seismic signals have affected the noise characteristic of GPS-derived displacement waveforms. S-transform is employed to assess the GPS capability to detect the seismic arrival time. The SM station AMT and the GPS station AMAT are in good agreement in seismic wave picking, and the difference is only 1.2?s in the north component, suggesting that the outcome of seismic wave picking using GPS data is reliable. Then, a classic empirical formula is employed to determine the moment magnitude. A robust moment magnitude (Mw 5.90) can be estimated by the nine GPS stations with about 23.9?s. If four GPS stations near the epicenter is chosen to determine the magnitude, it only take 13.0?s to retrieve a reliable preliminary (Mw 5.82) magnitude, which is 5.4?s ahead of nine stations. In addition, Cross Wavelet Transform (XWT) is adopted to measuring the correlation and phase relationship between GPS and SM records. The result of XWT analysis indicates 10?Hz GPS is capable of capturing reliable and accurate coseismic dynamic deformations, as evidenced by the XWT-based semblance being close to 1 between GPS and SM records. The above results confirm the capability of 10?Hz GPS to capture coseismic dynamic deformations, detect seismic arrival time, and determine earthquake magnitude. Moreover, rapid magnitude determination based on 10?Hz GPS data can be regarded as an important supplement to Earthquake Early Warning (EEW).